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Flatness-based Deformation Control of an Euler-Bemoulli Beam with 

In-domain Actuation 
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Abstract 

This paper addresses the problem of deformation control of an Euler-Bemoulli beam with in-domain actuation. The proposed 
control scheme consists in first relating the system model described by an inhomogeneous partial differential equation to a target 
system under a standard boundary control form. Then, a combination of closed-loop feedback control and flatness-based motion 
planning is used for stabilizing the closed-loop system around reference trajectories. The validity of the proposed method is 
assessed through well-posedness and stability analysis of the considered systems. The performance of the developed control 
scheme is demonstrated through numerical simulations of a representative micro-beam. 


I. Introduction 

The present work addresses the control of an Euler-Bernoulli beam with in-domain actuation described by an inhomogeneous 
partial differential equation (PDE). A motivating example of such a problem arises from the deformation control of micro¬ 
beams, shown in Eig.[T] This device is a simplified case of deformable micro-mirrors that are extensively used in adaptive optics 
El, El- Due to technological restrictions in the design, the fabrication and the operation of micro-devices, design methods that 
may lead to control structures with a large number of sensors and actuators are not applicable to microsystems with currently 
available technologies. 

One of the standard methods to deal with the control of inhomogeneous PDEs is to discretize the PDE model in space to 
obtain a system of lumped ordinary differential equations (ODEs) ID, E). Then, a variety of techniques developed for the 
control of finite-dimensional systems can be applied. However, in addition to the possible instability due to the phenomenon 
of spillover 0, the increase of modeling accuracy may lead to high-dimensional and complex feedback control structures, 
requiring a considerable number of actuators and sensors for the implementation. Therefore, it is of great interest to directly 
deal with the control of PDE models, which may result in control schemes with simple structures. 

There exists an extensive literature on the control of flexible beams described by PDEs in which the majority of the reported 
work deal with the stabilization problem by means of boundary control (see, e.g., 0, 0, 0, 0, 0, EHl). Placing actuators 
in the domain of the system will lead to inhomogeneous PDEs 0, ED, ED- The stabilization problem of Euler-Bernoulli 
beams with in-domain actuation is considered in, e.g., ED, E3, ED, ED, ESI- The deformation control of beams has gained 
in recent years an increasing attention and impotence (see, e.g., ED, ED, ED). 

In the present work, we employ the method of flatness-based control, which has been applied to a variety of infinite- 
dimensional systems (see, e.g., El, ED, ED, ED, El, El, El, Ei, ED). Nevertheless, applying this tool to systems 
controlled by multiple in-domain actuators leads essentially to a multiple-input multiple-output (MIMO) problem, which still 
remains a challenging topic. It is noticed that a scheme proposed recently in El, ED tackles this issue by utilizing the 
Weierstrass-factorized representation of the spectrum of the input-output dynamics. Nevertheless, an early truncation is still 
required in order to obtain a finite-dimensional input-output map, which is essential to apply the technique of flat systems. 

The control scheme developed in this paper is based on an approach aimed at avoiding early truncations in control 
design procedure. Specifically, it consists in first relating the original inhomogeneous model with pointwise actuation to a 
system in a standard boundary control form to which the technique of flatness-based motion planning may be applied for 
feedforward control. By using the technique of lifting, this boundary-controlled system can be transformed into a regularized 
inhomogeneous system driven by sufficiently smooth functions. It is shown that the steady-state solution of the regularized 
system can approximate that of the original system. A standard closed-loop feedback control is used to stabilize the original 
inhomogeneous system. Moreover, the MIMO issue is solved by splitting the reference trajectory into a set of sub-trajectories 
based on an essential property of the Green’s function. Well-posedness and stability analysis of the considered systems is also 
carried out, which is essential for the validation of the developed control scheme. It is worth noting that the approach proposed 
in this paper does not involve any early truncation. Moreover, it is shown that the invertibility of the resulting input-output 
map can be guaranteed a priori. This allows for the control design to be carried out directly with the original PDE model in 
a systematic manner. 
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Fig. 1. Schematic of the deformable microbeam. 


The remainder of the paper is organized as follows. Section describes the model of a representative in-domain actuated 
deformable beam. Section [III] introduces the relationship between the interior actuation and the boundary control and examines 
the well-posedness of the considered systems. Section [TV] addresses the feedback control and stability issues. Section [V] deals 
with motion planning for feedforward control. A simulation study is carried out in Section VI and, finally, some concluding 


remarks are presented in Section VII 


II. Model of the In-domain Actuated Euler-Bernoulli Beam 

As shown in Fig. [T] the considered structure consists of a continuous flexible beam and an array of actuators connected to 
the beam via rigid spots. As the dimension of the spots connecting the actuators to the continuous surface is much smaller 
than the extent of the beam, the effect of the force generated by actuators can be considered as a pointwise control represented 
by Dirac delta functions concentrating at rigid spots. Furthermore, it is supposed that the structure is supported at the left 
end, terminated to the frame, and is suspended at the right end, vertically suspended on a moving actuator, corresponding to 
a simply supported-shear hinged configuration HE). 

As the developed scheme uses one actuator for stabilizing feedback control, for notational convenience we consider a setup 
with N + 1 actuators. The displacement of the beam at the position x and the instance t is denoted by w{x, t). The derivatives 
of w with respect to its variables are denoted by and Wt, respectively. Suppose that the mass density and the flexural 
rigidity of the beam are constant. Then, the dynamic transversal displacement of the beam with pointwise actuators located at 
{xi,X2, ■ ■ ■ ,xn,xn+i} in a normalized coordinate can be described by the following PDF iflTll . ifTSll . ESll : 

N+l 

Wu{x,t) + Wa;xxx{x,t) = ^ aj{t)5{x - Xj), x G (0,1) ,t > 0, (la) 

i=i 

w{0,t) = = 0, f > 0, (lb) 

w{x, 0) = hQ{x),Wt{x, 0) = hi{x), X G (0,1), (Ic) 


where x is a normalized variable spanned over the domain (0,1), 5{x — Xj) is the Dirac mass concentrated at the point 
Xj G (0,1] and aj : t M., j = 1,N + 1, are the control signals. Hq and hi represent the initial value of the beam. 
Without loss of generality, we assume that 0 < xi < X 2 < • • • < xn < xn+i < 1- 

Remark 1: The considered problem can also be modeled as a serially connected beam ca. Indeed, these two models are 
equivalent in the sense that they lead to the same abstract linear system in variational form (see, e.g., ESI). 

The control objective is to steer the beam to a desired form in steady state via in-domain pointwise actuation. 

The model given in ([T]) is a nonstandard PDE due to the unbounded inputs, represented by Dirac functions, on the right-hand 
side of ( [Ta| . Therefore, we need to invoke the weak derivative in the theory of distributions 130) and discuss solutions of the 
considered problem in a weak sense. To this end, we start by defining the convergence on C^{n). 

Definition 1: Fet fl be a domain in K. A sequence {pj} of functions belonging to C(j°(n) converge to tp € CQ°(n), if 

(i) there exists K C D. such that supp((/jj — ip) C K for every j, and 

(ii) limj_>oo 9"v^j(x) = d'^(p{x) uniformly on K for all n > 0. 

The linear space having the above property of convergence is called fundamental space, denoted by The 

space of all linear continuous functionals on denoted by ^'(17), is called the space of (Schwartz) distributions on 17, 

which is the dual of ^(17) (see, e.g.. Chapter 1 of lISTl for more properties of ^(17) and ^'(17)). 

Denote the set $ by 

$ = {</> e iT2(0,l);(/.(0) = </>,(!) =0}. (2) 

We give the definition of weak solution of Q- 

Definition 2: Fet T > 0 and ai G (0,T) for i = 1,2, ...,7V + 1. Fet Tig € hi G (0,1). A weak solution to the 
problem Q is a function w € (^([OjT]; $) H C^([0,T] ;L^(0,1)), satisfying 

w{x, 0) = ho{x), X G (0,1), 
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such that, for every v G <^^([0, T] ; ^(0,1)), one has for almost every t S [0, T] 

/ Wt{x,T)v{x,T)dx — / hi{x)v{x,0)dx 

Jo Jo 


oT /*1 


Wt{x^ t)vt{xj t)dxdt 


/o ^0 

pr rl 


Wxx{x,t)vxx{x,t)dxdt 


JO Jo 

N-\-l rl 


= E 


ai{t)S(x — Xi)v(x, f)dxdf. 


/o ^0 


(3) 


III. In-domain Actuation Design via Boundary Control 
A. Relating In-domain Actuation to Boundary Control 

In the scheme developed in the present work, we use N feedforward control signals to deform the beam, while dedicating 
the feedback stabilizing control to one actuator. For notational simplicity, we assign the feedback control to the {N + 1)* 
actuator. Hence, the dynamic of the desired trajectory driven by the feedforward control can be expressed as 

N 

wft{x,t)+wi^^^{x,t)=J2<^ji0Hx-Xj), xG{0,l),t>0, (4a) 

= 0, t>0, (4b) 

w‘^{x, 0) = wf{x, 0) = 0, x G (0,1). (4c) 

A weak solution of 0 can be defined in a similar way as that of Q described in Definition Note that it is shown later that 
the initial condition given in © will be captured by the regulation error dynamics. Therefore, feedforward control design can 
be carried out based on System Q with zero-initial conditions. 

Due to the fact that the model given in (j^ is driven by unbounded inputs, we will apply a sequence of cutoff functions 
(p{x — Xj) having unit mass and satisfying (x — Xj)(p(x — Xj)dx = 0, also called blobs Il32]| . to approximate 6{x — Xj) in 
the sense of distributions. We will then show that in steady state, w‘^ can be approximated by a sufficiently smooth function. 
Specifically, we consider first the following boundary controlled PDF; 


Utt{x,t)u^xxx{x,t) = 0, xG{0,l),t>0, (5a) 

u{0,t) = = 0, t>0, (5b) 

UxxxO,t) = 9(0^ (5c) 

u(a:, 0) = ut(a;, 0) = 0, x S (0,1), (5d) 


where g(t) = 9ji0- Throughout this paper, without special statements, we assume that gj{t) G C^([0,-|-c»)) and 

Pj(0) = gj{0) = 0 for j = 1,2, Notice that the motivation behind considering the system of the form given in (|^ is 

that it allows employing the techniques of boundary control for feedforward control design while avoiding early truncations 
of dynamic model and/or controller structure. 

A weak solution to (|^ can be defined in a similar way as that of ([T]l described in Definition S 

Definition 3: Let T > 0. A weak solution to the problem 0 is a function u G Cjp, T]; $)n^([0, T] ; T^(0,1)) satisfying 

m(x, 0) = itt(x, 0) = 0, x G (0,1), 


such that, for every v G C^([0,T] ;<!)), one has for almost every t G [0,T] 


/ Ut(x, t)v(x, T)dx — / / Ut(x,t)vt(x,t)dxdt 

Jo Jo Jo 


g{t)v{l,t)dt-\- / Ua:x{x,t)Vxa;{x,t)dxdt = 0. 

Jo Jo 

Let 'fi{x, t) = u{x, t) — 9j{0Hj{x), where Hj{x), j = 1, 2,..., N, are defined in [0,1] satisfying 



Hjxxxx{x) — ^{x Xjfi 

iT,(0) = = 0, = 1, 


(6) 


(7a) 

(7b) 
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with </j(a;) G L^(IR). Then by lifting, Q can be transformed into the following PDE with zero boundary conditions 


iptt{x,t) + V' xxxx (x,t) 

N N 

= - '^gj{t)Hj{x) - a; G (0,1),f > 0, (8a) 

^/>(0,f) = = il^xxi0,t) = ipxxxi^,t) = 0, (8b) 

^/>(a;, 0) ='0t(x,O) = 0, x G (0,1). (8c) 

In order to relate 0 to 0, we hrst establish a relationship between Q and especially in steady state. Let aj{t) = —gj{t) 


and suppose that limt_>.oo ctj{t) = aj and limt_).oo gj{t) = g^ for j = 1,..., N. We have then in steady state; 

N 

i^xxxxix) = ~ ^ '^ g-jHjxxxx{x) 

N 

= ^ ajip{x — Xj), X € (0,1), f > 0, (9a) 

i=i 

V^(0) = ^,(1) = ^,,(0) = ^_(1) = 0. (9b) 


For j = 1,2, N, given a sequence of blobs {(pmix — Xj)}, we seek a sequence of functions {Hp{x)} such that 


H^xxx (X) = (f^{x - Xj) (10) 

with H'p{x) satisfying and (pm{x — Xj) 6{x — Xj) in ^'(0,1) as m —+oo. Therefore, considering 0 and the 
steady-state model of 0. we have 'ijj^ —> w'^ in as m —> +(x; (see Theorem]^. Hence, for feedforward control 

design, we may consider the systems ([^ and 0. 

Lemma 1: f32\ Let (pm(x) G L^(IR) be dehned by 


(Pmix) 


1 sinmx 

TT X 


( 11 ) 


Then ipm{x) has the following properties: 

, , f~^°° 1 sinma: 

(i) / ipm[x)dx = /-dx = 1. 

J — oo J — oo ^ ^ 

(ii) ipm{x) —>■ 5{x) in ^'(M) as m —>■ oo. 

By taking Lpm given in ([TT]l as the input to ([TOll, we get 





Xj)dsdtdydz. 


Theorem 2: Let (fim{x) be dehned as in Lemma Assume that aj{t) = —gj{t) tends to aj = —gj as f —>■ oo for all 
j = 1,2, Denote by tpj and by Wj the steady-state solutions of System Q and System ([^, respectively. Then tpj Wj 

in 1]) as m ^ -foo for j = 1, 2,..., N. 

Proof: In the steady state, we have 

jxxxx^^^ 9j ^jxxxxi^) CtjPmix Xj), (12a) 

^7(0) = (1) = = 0, (12b) 

and 


Taking v{x) G S’{0, 1), with ti(0) 


Wjxxxxix) = ajSix-Xj), 

w^(0) = u;7(l) = wf,^x(0) = w^xxxW = 0- 
'fa;(l) = XxxiO) = 0, as a test function and integrating by parts, we get 



— 



X — Xj) — S(x 


Xj)) v(x)dx. 


(13a) 

(13b) 
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Since ipm{x — Xj) —>■ 6{x — Xj) in the sense of distributions as m ^ +oo and Vxxx G ^(0,1), it follows that —>■ 

in the sense of distributions as m —)■ +oo for j = 1, 2, N. Furthermore, as G L^(0,1) and G L^(0,1), we have 
1 /)^ -G- a.e. in (0,1) (see Lemma 3.31 of OTI . page 74). Then by the continuity of i/>j^ and Wj^, we conclude that 
—)■ Wjx pointwise in (0,1). Therefore —>• Wj in C'^([0,1]) as m —>■ +oo for j = 1, 2,..., N. ■ 

Remark 2: Note that describing the control acting on discrete points in the domain by Dirac delta functions is a mathematical 
abstraction, which allows for the beam equation to possess some fundamental properties, in particular the exact controllability 
(see, e.g., |l33l), and facilitates control design. However, the real physical devices cannot generate unbounded control concen¬ 
trated on a point. From this viewpoint. Theorem provides an assessment of the gap between the performance achievable by 
utilizing real physical devices and that predicted by the corresponding idealized mathematical model. 


B. Well-posedness of Cauchy Problems 

Well-posedness analysis is essential to the approach developed in this work. In this subsection, we establish the existence 
and the uniqueness of weak solutions of equations Q, 0, 0 and 0. 

Theorem 3: The following statements hold true; 

(i) Assume aj G L^(0,T) for j = 1,2, ...,N + 1. Let hg S G L^(0,1) and T > 0. Then System 0 and System 

0 has a unique weak solution w G C'([0, T]; $) n T] ; L^(0, 1)) and w'^ G T]; $) n C'^([0, T] ; L^(0, 1)), 

respectively. 

(ii) Let T > 0 and g G C'^([0,r]). Then System 0 has a unique weak solution u G (7([0,T]; T*) n C'^([0,T] ;L^(0, 1)). 
Furthermore, if p G C^([0, T]) and H'^{x), j = 1, 2,..., N, is defined as in ( |7bl l, then System 0 has a unique solution 
V'™eC([0,T];$)nCi([0,T];L2(o,l)). 

Proof: The proof of (i) can be proceeded step by step as in Proposition 3.1 of im. We prove the first result of (ii) and 
the second part can be proceeded in the same way. We assume first g{f) G (^^([OjT]) and consider the following system: 

vtt{x,t) + Vxxxx{x,t) = Qa; - g{t), x G n,t > 0, (14a) 

1/(0,f) = Vxii,t) = i/xx(0,f) = iyxxx{i,t) = 0, (14b) 

iy{x,0) = 0,ot{x,0) = 0, X G Cl. (14c) 


Let X = $ X L^(0,1) and iJpjj(0,1) = {u G H‘^{0,l)',u{0) = Ux{l) = Uxx{0) = Uxxx{^) = 0}. Define the inner product 

on X by ((ui,ui), (^ 2 ,W 2 ))x = fgiuixxUsxx + viV 2 )dx. Define the subspace T>(Ao) C X by T>(Ao) = {(u,v); (u,v) G 
X 4)}, with the corresponding operator Ag: ^{Ag) -G X defined as: 


/4o 




One may easily check that V^Ag) is dense in X, Ag is closed, Ag = —Ag, and {Agz, z)x = 0. Thus, by Stone’s theorem, Ag 
generates a semigroup of isometries on X. Based on a classical result on perturbations of linear evolution equations (see, e.g.. 
Theorem 1.5, Chapter 6, page 187, ll34ll ) there exists uniquely z = (zi,Z 2 ) G (^^([OjT] ;Ar) nC([0,T] ;'D{Ag)), such that 


dz 

dt 



z(.,0) = (0, 0), 


which implies that ( [T4l l has a unique solution v = Zi G C{[0,T] 1)) n C'^([0,T]; $) in the usual sense. Particularly, 

J/ G C([0, T]; 4>) n C^([0, T ]; L^(0,1)) is a weak solution. A direct computation shows that u = v — (^a; — g{t) is a 
solution of 0 in the usual sense and, in particular, it is a weak solution. 

Now for g G C'^([0,r]), let S C^([0,T]) such that —>• p in C'^([0,T]). Consider = i/„ — (^x — ^x^) p„(f), where 
i/n is the solution of ( [T4| ) corresponding to the data gn{t). Then arguing as above and taking limit as in Chapter 2 of llT5l . we 
can obtain the existence and the uniqueness of a weak solution of 0. ■ 


IV. Feedback Control and Stability of the Inhomogeneous System 

The validity of the proposed scheme requires a suitable closed-loop control that can guarantee the stability of the original 
inhomogeneous system. As the {N + 1)* actuator is dedicated to stabilizing control, ( [Ta| l can be written as 

Wtt{x,t) + w 

XXXX {x, t) - aN+i{t)5{x - XN+i) 

N 

= aj{t)S{x — Xj), X G (0,1), f > 0. 


(15) 
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Suppose further that the feedback control is taken as im, M 


aN+i{t) =-kwt{xN+i,t), (16) 

where fc is a positive-valued constant. Then in closed-loop, becomes 

Wtt{x,t) + w 

XXXX {x, t) + kwt (x, t) S(x — xn+i) 

N 

= aj{t)d{x — Xj), X £ {0,1), t > 0. (17) 

i=i 

Let 'D{A) = {{w, v);{w,v) G {H^{0, l)n (7L‘(0, a;i)U77^(a;i, a; 2 )U • • ■UH'^{xn, xn+i)^H'^{xn+i, 1))) x 77^(0,1), ti;(0) = 
Wx{l) = Wxx{0) = 'Wxxx{l) = 0,'(;(0) = Ua;(l) = 0} and X be dehned as in the proof of Theorem]^ Let T > 0 and 
Oij G L^(0,T) for j = 1,...,X. Assume {ho,hi) G X. To address the stability of System ( [TT] ) with the boundary conditions 
and initial conditions we consider the corresponding linear control system under the following abstract form: 


i = Az + Ba, t > 0, 

z{0) = z° = {ho,hi)'^, 

where z = {w, v)'^, a = (ai, • • • , A: I?(A) —>■ X is dehned as: 


A 


\ '^xxxx kv6{x 

with V = Wt and B: -G 'D'{A*), where A* is the adjoint of A, which is dehned as: 

“ ' y, \JyGV{A*). 




Note that B*\ V{A*) 


pN 


is dehned by 


(18a) 

(18b) 


(19) 


( 20 ) 


B*y={y2{xi),...,y2{xN))'^,yy= GV{A*). (21) 

Let U = . Then the solution of the Cauchy problem ( [T8] l can be dehned as follows (see llT5l . Dehnition 2.36, p53): 

{z{T),y'^)x 

={z°,S*{T)y^)x+ na{t),B*S*{T-t)y^)udt, (22) 

Jo 

y TG[0,T],y y^ GX. 

One may verify, as in Chapter 2 of 051 . that the solution dehned by ( |2^ is also a weak solution of under the form given 
in Dehnition 12] 

Theorem 4: Assume aj G L°°(0,-t-oo), j = l,...,iV. For any xx+i G (0,1] being a rational number with coprime 
factorization, there exist positive constants Ci, C 2 and A such that for any 0 < r < T < 00 , there holds: 


||^(t)||x < Cie-^nk°llx + C2||a|U=o(0,T). (23) 

Proof: The proof is a slight modihcation of Theorem 2.37 in 051 . First, the admissible property of B can be obtained as 
(3.38) in HU. A direct computation gives 

)x “ 

Therefore, A is a dissipative operator. Moreover, it can be directly verihed that A is onto and hence, according to Theorems 
4.3 and 4.6 of EQI (p. 14-15), it generates a Co-semigroup of linear contractions S(t) acting on X. Furthermore, S(t) is 
exponentially stable if xx+i G (0,1) is a rational number with coprime factorization, in particular for xx+i = 1 (see, e.g., 
CD, US), i.e. there exist two positive constants Ci and A, such that 


ll>5'(f)||£(JV;X) < C'le V t > 0. 
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Then for a weak solution defined in ( [22| ), we have that for all 0 < r < T < cx): 

{z{T),y^)x 

= {S{T)z°,y^)x + f {S{T-t)Ba{t),y^)udt 
Jo 

= {S{T)z°,y^)x + l^J S{T-t)Ba{t)dt,y^'^ 

<Cie-"nk°llx||2/1|x + Ci||a|Uo.(o,T)b1|L»(o.r) 

Jo 

<C,e-^n\z°\\x\\yl\x + ^||a|U~(o,T)||?/1|L<» 

<C,e-^n\z°\\x\\yl\x + ^\\ah<^^o,T)\\yl\m- 

Applying the interpolation formula 

||tta:||L 2 ^ (l|tt||L 2 “t“ IIl^) 5 

where c is a positive constant, and the Sobolev embedding theorem yields 

llVIlrri < c||Vllr/2 < c\\y^\\x- 

We have then 

{z{r),ynx < {c,e-^^\\z^\\x + ^||a|U=o(o,T)j IIVIU, 

which implies that ( |23l l holds with = CicjX. ■ 

Now consider the regulation error defined as e{x, t) = w{x, t)—w‘^{x). Denoting Aaj{x, t) = aj{x, t)—aj{x), j = 1,. .., N, 
then from ( [TtI i, ( [Tbl l, ( [Tc| l and the steady-state model of Q, the regulation error dynamics satisfy 

ett{x,t) + 

^XXXX {x, t) -I- ket {x, t) S(x — xx+i) 

N 

= Aaj{t)S{x — Xj), X € (0,1), t > 0, (24a) 

i=i 

e{0,t) = =0, t > 0, (24b) 

e(a;, 0) = eo{x) = ho{x) — w^^^x), 

et{x,0) = ei{x) = hi{x), xG( 0, 1). (24c) 

Obviously, the regulation error dynamics are in an identical form as ( [TtI i with the same type of boundary conditions. We 
can then consider the solution of System ( |24| ) defined in the same form given by ( |22l l with Zg = (e, et) and z® = (eo,ei). 

Corollary 5: Assume that all the conditions in Theoremare fulfilled and z® G X. Then there exist positive constants Ci, 
C 2 and A, independent of t, such that for any f > 0, there holds: 

||^e(t)||x < Cie-^iz^llx + C2||Aa|U=o(o,+oo)- (25) 

Moreover, if limt_>oo Aaj = 0 for all j = 1,..., N, then limt_).oo e{x, t) = 0, Va; G (0,1). 

The feedforward control satisfying the conditions for closed-loop stability and regulation error convergence can be obtained 
through motion planning, as presented in the next section. 

V. Motion Planning and Feedforward Control 

According to the principle of superposition for linear systems, we consider in feedforward control design the dynamics of 
System 0 corresponding to the input gj{t) described in the following form: 


“ 1 “ '^jxxxxis^ — 0 ; 

(26a) 


(26b) 

(I 7 

(26c) 

Uj (x, 0) = Ujt (x, 0) = 0. 

(26d) 


The required control signal gjit) should be designed so that the output of System ( |2^ follows a prescribed function Uj{xj,t), 
while guaranteeing that the conditions in Theorem and Corollary are fulfilled. 


Motion planning amounts then to taking Uj{xj,t) as the desired output and to generating the full-state trajectory of the 


subsystem Uj{x,t). By applying a standard Laplace transform-based procedure (see, e.g., ITSl . ||22l, |l27l), we can obtain 
the full-state trajectory with zero initial values expressed in terms of the so-called flat output, yj{t), and its time-derivatives: 


CO / n 

+ E E 




n—l \k—Q 


(4fc -f l)!(4(n—/c) -I- 2)! 


(27) 


-E 


j,4fc+3 


(4fc-f 3)!(4(n-fc))! 




Now let yj{t) = yj(j>j{t), where is a smooth function evolving from 0 to 1. The full-state trajectory can then be written 


. N - (1 Is 

Uj(x,t) =y ■ [-X- -X 


iit) 


2 6 
oo / n 

E (^E (4A: + l)!(4(n-A:) + 2)! 


X 


Ak+l 


(28) 


-E 


X 


,4/c+3 


(4A: + 3)!(4(n-A:))! 




The corresponding input can be computed from P6c| l, which yields 

OO / n 

9j{t) = - E (E 


n—l \k—l 


-E 

k=0 


(4fc)!(4(n-fc))! 


(4fc- 2)!(4(n-fc) -f 2)! 


(29) 


For set-point control, we need an appropriate class of trajectories enabling a rest-to-rest evolution of the system. A convenient 
choice of (pj{t) is the following smooth function: 


0 , 


= < 


[ exp(-l/(T(l - 

Jo _ 

[ exp(-l/(T(l - rjjj'^dr 

Jo 


1 , 


if f < 0 


if f G (o,r) 


if f > T 


(30) 


which is known as Gevrey function of order a = \ + l/e, e>Q (see, e.g., ||27l). 

For the convergence of ( |28| ) and ( |29| ), we have 

Proposition 6: If 4>j{t) in the basic output y^{t) = y^(j)j{t) is chosen as a Gevrey function of order 1 < ct < 2, then the 
infinite series ( |28] l and ( |29] l are convergent. 

The proof of Proposition is based on the bounds of Gevrey functions of order a given by (see, e.g., 121], 1^ ): 

{kir 


3K,M> 0,Vfc G Z>o,Vf G [to,T ], 


< M- 


Rk ■ 


(31) 


The detailed convergence analysis is given in Appendix |IX-A| 

Remark 3: In general, the Gevrey bounds are unknown, but can be estimated following the way presented in 11211 . Further¬ 
more, a symmetric function in the transient phase can be considered to improve convergency analysis M- 
Now let 

, , 1 1 o 

Pj{x) =2^~ 0 ^ ’ 

fX nZ py nt 


nX nZ ny nl 

Ij,m{x) = / / / PmO - Xj)dsdtdydz, 

Jo Ji Jo Ji 

^j.n(a;) =(^E (4fc+ l)!(4(n-fc) -f 2)! “ (4fc-f 3)!(4(n-A:))! 




1 


(4fc - 2)!(4(n-fc) -f 2)! ^ (4fc)!(4(n-A:))! 


-E 


1 


(-ir, 

(-ir- 
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Set ip'^{x,t) = iprn{x,t),m > 0. By the definition of ij){x,t), we obtain 




n—1 


( 33 ) 


To compute Ij^rn{x), we note that 


If, 


,{x-x,) = - ^ 


(- 1 ) 


fc+1^2fe —1 


(2fc-l)! 


-(X-Xj) 


2k-2 


Therefore, 


^ oo 

=— 'y' 

TT • * 


(_l)fc+1^2fe-l ^( 2 ; _ 

(2fc + 2)! 


^ (2fc - 1) 


(1 - Xjf 


xj^x"^ (1 - 


6((2/c-l)!) 2(2fc)! (2fc + l)! 


(34) 


x'^’^x (1 - Xj) 


+ 


2fc-l„ 


„2/c+2 


(2fc)! 2(2fc-1)! (2fc + 2)!y' 


Claim 7: given in ([3^ is convergent for all x,Xj G (0,1) with respect to any fixed m. 

The proof of this claim is given in Appendix |IX-B| 

Based on Claim the series in the right hand side of are convergent and tjjJ^{x,t) can be expanded by and ( [34l l. 
Moreover, for gj{t) given in ( |29l l, ij}J^{x,t) tends to (x) — y^Ij^mix) as f —>■ 00 . Note that for fixed x, the radius of 
convergence of ipj ^ on m is 00 . Thus, we can let m —?> + 00 . 

To complete the control design, we need to determine the amplitude of flat outputs yj, j = 1,... ,N, from the desired shape 
w‘^{x) that may not necessarily be a solution of the steady-state beam equation ( pB] ). To this end, we propose to interpolate 
w^{x) using the Green’s function of 0, denoted by G(x,^), which is the impulse response corresponding to the input 
5{x — ^) and is given by 

r T-a / \ 

0 < X < 




6 


j < X < 1. 


Due to the principle of superposition for linear systems, the solution to ([T3|l, w^{x), can be expressed 


as 


pi N 

w‘^{x) = Y 
Jo .-I 


^0 

N 


= YG{x,Qaj. 

3=1 

Now taking N points on w‘^{x) and letting W'^{xj) = w‘^{xi = j), j = 1,..., N, yield 

'G(xi,Ci) ... G(xAr,^i)\ /“i\ fw'^{xi)' 

^G(xi,^Ar) ... G{xn,^n)/ \aN/ \w'^{x]v)^ 

which represents a steady-state input to output map. 

Claim 8: The map of ( |T5| l is invertible for all Xj,^j G (0, 1), j = 1,..., N and Xi ^ Xj,^i ^ if i ^ j. 
The proof of this claim is given in Appendix IX-C| 

As aj = —gj and limt_>.oo gj(t) = gj = —y^ for all j = 1,..., N, we obtain from Claim that 


(35) 


Vi 


kVn/ 


/G(xi,^i) ... G{xN,ii)\ w‘^{xi) 


\G(xi,^Ar) ... G{xn,^n)/ \w'^{xn)^ 


( 36 ) 
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VI. Simulation Study 

In the simulation study, we consider the deformation control in which the desired shape is given by 

w\x) = - 10-3 

^3^_4oo(x-o. 7)=^ , cce (0,1), 


(37) 


as shown in Fig. 2(a) Note that the maximum amplitude of the shape given by ( |J7] ) is 3.8 x lO-^, which represents a typical 
micro-structure for which the beam length is of few centimeters and the displacement is of micrometer order. 

In order to obtain an exponential closed-loop convergency, the actuator for feedback stabilization is located at the position 
xn +1 = 1. To evaluate the effect of the number of actuators to interpolation accuracy, measured by — ui'^(a:)||Li(o,i) 

and control effort, we considered 3 setups with, respectively, 8, 12 and 16 actuators evenly distributed in the domain. It can 
be seen from Fig. that the setup with 8 actuators exhibits an important interpolation error and the one with 16 actuators 
requires a high control effort in spite of a high interpolation accuracy. The setup with 12 actuators provides an appropriate 
trade-off between the interpolation accuracy and the required control effort, which is used in control algorithm validation. 

A MATLAB Toolbox for dynamic Euler-Bernoulli beams simulation provided in Chapter 14 of llJTll is used in numerical 
implementation. With this Toolbox, the simulation accuracy can be adjusted by choosing the number of modes used in 
implementation. In the simulation, we implement System (lljl with initial conditions ho(x, 0) = —3 x I0“3g-40o(x-o.8) 
and hi{x,Q) = 0. As an undamped beam is unstable in open-loop, the controller tuning is started by determining a suitable 
value of the closed-loop control gain k. The basic outputs (pj{t) used in the simulation are Gevrey functions of the same 
order. To meet the convergence condition given in Proposition]^ the parameter in ( |30| is set to e = 1.111. The corresponding 
feedforward control signals with T = 5, ai,..., ai 2 , that steer the beam to deform are illustrated in Fig. The evolution of 
beam shapes and the regulation error are depicted in Fig. It can be seen that the beam is deformed to the desired shape 
and the regulation error vanishes along the entire beam, which confirms the expected performance of the developed control 
scheme. 




(a) (b) 

Fig. 2. Effect of number of actuators: (a) interpolation accuracy {w‘^{x): desired shape; w'^{x): solution of the steady-state beam equation with n =8, 12 
and 16 in-domain actuators); (b) amplitude of steady-state control signals for different setups. 



Fig. 3. Feedforward control signals: ai,... ,ai 2 . 


VII. Conclusion 

This paper presented a solution to the problem of in-domain control of a deformable beam described by an inhomogeneous 
PDF. A relationship between the original model and a system expressed in a standard boundary control form has been 
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established. Flatness-based motion planning and feedforward control are then employed to explore the degrees-of-freedom 
offered by the system, while a closed-loop control is used to stabilize the system around the reference trajectories. The validity 
of the developed approach has been assessed through well-posedness and stability analysis. System performance is evaluated by 
numerical simulations, which conhrm the applicability of the proposed approach. Note that in order to improve the resolution 
of manipulation, more actuators may be expected. Nevertheless, for the control of one-dimensional beams, the proposed scheme 
contains only one closed loop, allowing a drastic simplihcation of control system implementation and operation. This is an 
important feature for practical applications, such as the control of large-scale deformable micro-mirrors in adaptive optics 
systems. 
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IX. Appendices 

A. Proof of Proposition 

We prove the convergence of the power series ( |28| l using Cauchy-Hadamard Theorem. The convergence of 
easily using the same argument. 

Denote in ( |27| ): 


then follows 


bn 


„4fc+l 


Kk=0 


(4fc + l)!(4(n—A:) + 2)! 


-E 

fc=0 


j,4/c-|-3 


(4fc + 3)!(4(n-fc))! 


i-irfYit). 


Then, ( |28] l converges if limsup„_^oo ^\bn\ < 1- 

Thus, for 0 < Xj < 1 and x G [0,1], can be majorized as: 


\bn\ < 
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E 

n 

E 

k^O 

n 

E 

n 

E 


1 


(4fc + l)!(4(n-fc) + 2)! 
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(4A + 3)!(4(n-A:))! 
4fc + 2 




(4fc + 2)!(4(n-fc))! 


4fc + 2 


k=0 

<2(4n + 2) 


(4A: + 2)!(4(n-A:))! ' 
1 




<2(4n + 2)- 


E 

k=0 

16 


(4A + 2)!(4(n-A:))! 


fC\t) 


Therefore, 


(4n + 2)! 

lim sup < lim sup 16 


^Y\t) 


< 


2M 16”((2n)!)'" 
(4n+l)! 

^2My/" ((2n)!)'"/” 


= —z lim sup 


yx'^nj ((4n + 1)!)V’T’ 

(((2n)!)i/(2»))2Wn 


K'^ n-s-oo (((4n+l)!)i/(4"+i))(‘^"+i)/"' 


Applying Stirling’s formula s/nJ ~ (n/e) yields: 

lim sup < 






(38) 


( 39 ) 
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Since 


{2nfr,alr. 


‘“-.r (4n+ 1)<1"W/1 


= lim sup 


(2n) 


n—>00 (4ri. + 



(40) 


we can conclude by Cauchy-Hadamard Theorem that ( [28] l converges for a < 2 with a radius of convergence of infinity and 
for CT = 2 if > 1. The series (|28ll diverges if ct > 2. 


B. Proof of Claim 

Consider the first series in ([34ll. Fixing m > 0, for x, Xj G (0,1), j = 1, 2,..., TV, we have 


(_l)fe+1^2fc 1 ( 2 ; _ 2.^.^2fc+2 


( 2 fc-l) ( 2 fc + 2 )! 
Since lim„_>.oo ^=7 = 0, it follows that for any a > 0 


< 


,2fe-l 


(2fc- 1)!' 


/a" 1 

lim \ —r — ^ = = 0 . 

n—)-oo y 77,! n—>-oo ^ji\ 


Thus iTfc-i'i! convergent. We conclude that the first series in ( [34l l is uniformly convergent. The convergence of the 

other terms in ([3?| can be proved in the same way. 


C. Proof of Claim 

We denoted by CjOIatxat matrix in the left-hand-side of ( [T5] l formed by Green’s functions. We argue by contra¬ 

diction. If, otherwise, [G{xi,^j)]j^^^ is not invertible, then it is of rank less than N. Without loss of generality, assume that 
there exist TV — 1 constants, ki, k 2 , ■■■, kM-i, such that 

N-l 

G(a;i,eiv)= 


Z =1 

N-l 


G(a;2,etv)= ^fc4G(ir2,e4), 


Z =1 


The above equations show that 


N-l 

G{xn,^n) = ^ hG{xNi^i)- 

i^l 

A^-1 

G(x,C7v) = 


(41) 


has TV different positive solutions xi, X 2 , • ■ ■ Xn, Xi € {0,1), i = , N. 

We consider two cases: 

(i) If TV > 3, since T = 1,..., TV, are distinguished, G{x, ^j), i = 1,..., TV, are all different from each other. Hence 


G{x,^n) ^ fciG(a::,^i) -f k 2 G{x,^ 2 ) H - f kN-iG{x,^N-i)- 

Note that G{x, ^j),j = 1,..., TV, are of order at most 3, then ( |4T] l has at most 3 different solutions in M, which is a contradiction, 
(ii) If TV < 3, it is easy to check that ig has a solution x = 0, and a pair of solutions x = x^ and x = —x^ near the origin 
0. By the assumption, ( |4l] l has TV different positive solutions, then it must be TV = 1, which leads to a contradiction with the 
non-invertible property of G(a;i,^i) f 0 . 










